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Abstract 

We have analyzed the IIB matrix model on the basis of the improved mean 
held approximation (IMFA) and have obtained evidence suggesting that the four¬ 
dimensional space-time appears as its most stable vacuum. This method is a sys¬ 
tematic way to obtain an improved perturbation series and was first applied to the 
IIB matrix model by Nishimura and Sugino. In a previous paper, we reformulated 
this method and proposed a criterion for the convergence of the improved series, 
that is, the appearance of a “plateau.” In this paper, we carry out higher-order 
calculations, and find that our improved free energy tends to have a plateau, which 
shows that IMFA works well in the IIB matrix model. 
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1 Introduction and summary 


String theory has been proposed as a unified theory of fundamental interactions including 
gravity. If it is the well-defined quantum gravity theory which our nature adopts, it should 
be able to predict several properties of our universe, for instance, the gauge symmetry at low 
energy, the particle content and their masses, and the dimensionality of our space-time. The 
last of these is the property we would like to investigate in this paper. 


As is well known, superstring theory has infinitely many vacua with various dimensionalities 
that arise perturbatively. One of the most promising scenarios to construct a realistic four¬ 
dimensional model from superstring theory is a compactification in which the ten-dimensional 
space-time consists of a flat four-dimensional space-time and a six-dimensional compact man¬ 
ifold that is small enough to be invisible to our experiments. There are several ways to realize 
this scenario, for example, Calabi-Yau compactification ||, fermionic construction [[J, and so 
on, but they are all stable and there is no way to single out the true vacuum perturbatively. 
However, it has been revealed that this problem originates from perturbative formulations 
themselves, and to overcome this difficulty, we need a nonperturbative formulation and anal¬ 
ysis of string theory. 


In the mid-1990s, some models were presented as a constructive formulation of M-theory || 
and certain kinds of superstring theories, such as type IIA matrix string 0, type I superstring 

and heterotic superstring ||. Here we would like to analyze the model called the IKKT 
model ||, which is conjectured to be a nonperturbative definition of type IIB superstring 
theory. It seems to be most promising for nonperturbative analyses of superstring theory, and 
some kinds of extensions of this model have been proposed ]I(J. For a review of the IIB matrix 
model, see Ref. Also, there are some mechanisms proposed for dynamical breakdown of 
Lorentz symmetry in this model (see Ref. Hl6|). 

In order to analyse the IIB matrix model, we use a method that we call the “improved 
mean field approximation” (IMFA) [|T|. This approximation was applied to the large-A" reduced 
Yang-Mills models by Oda and Sugino [|f2|. Then, Nishimura and Sugino applied it to the IIB 
matrix model in an excellent work [|^]. They obtained a result that suggests the breakdown 
of Lorentz symmetry to the four-dimensional symmetry. In a previous paper ]l|, we analyzed 
how the IMFA scheme works and discovered a general structure that we call the “improved 
Taylor expansion.” Furthermore, we proposed a principle for choosing the mean fields, that 
is, the existence of a “plateau.” The emergence of a plateau indicates that the approximation 
scheme works well and, in fact, its existence can be confirmed in some exactly solvable models. 
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(See Ref. [l|] for reference. There, many examples are given which show how good this scheme 
is.) Furthermore we developed a computational method using two-particle irreducible (2PI) 
graphs, which simplifies the calculation drastically. The 2PI free energy has a close relationship 
with the Schwinger-Dyson equation as discussed in Ref. [|I|]. 

Using this method, we calculated the free energy up to 5th order and obtained a preliminary 
result which suggests that the eigenvalue distribution of the matrices preserves only the four¬ 
dimensional rotational symmetry. 

In this paper, we carry out a further calculation up to 7th order. We find evidence suggest¬ 
ing the emergence of a plateau, which was not clear at 5th order. In fact, we have evaluated 
the free energy and the extent of the eigenvalue distribution for various Lorentz symmetries by 
introducing corresponding mean fields. In the case of SO(4) symmetry, the number of extrema 
of the free energy increases as we go to higher orders, and it seems that they form a plateau. 
Furthermore, the eigenvalues are widely distributed in the four “non-compact” directions, 
while they tend to gather in the six compact directions. One the other hand, if we impose 
SO(7) symmetry, the number of extrema of the free energy does not grow enough to form a 
plateau, and the eigenvalues are distributed somewhat isotropically in the ten directions. As 
shown in Ref. [|lj, the other cases are reduced to the above two cases, or do not have a plateau 
at all. 

We evaluate the 2738 graphs to obtain the free energy up to 7th order. All the calculations, 
including the generation and computation of the graphs, are totally automatized now, and we 
will be able to go further and find reliable results. 

In §|2j. we provide a short review of the improved mean field approximation with a </> 4 matrix 
model as an example. We examine a distribution of extrema of the free energy of this model, 
which provides a method of the search for a plateau in the IIB matrix model. In we apply 
IMFA to the IIB matrix model. We obtain the free energy, investigate the distribution of its 
extrema to search for a plateau and examine the eigenvalue distributions. Section [| is devoted 
to a conclusion and discussion. In Appendix [A], we make use of a </> 4 -QED type matrix model 
to confirm our counting of the graphs. 
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2 Improved Mean Field Approximation for the Free En¬ 


ergy 


In this section, we review the improved mean field approximation that was developed in our 
previous paper |lj]. Although all the techniques we introduce here are the same as those used 
in Ref. [IJ, we explain them to make this paper self-contained and for reader’s convenience. 

2.1 Ordinary vs. Improved perturbation theory 

Suppose we have some action function S(x ) = S 0 (x ) + S i(x) and its free energy 



( 2 . 1 ) 


where S 0 (x) is the unperturbed part, which can be integrated analytically, and Si(x) is difficult 
or impossible to integrate and is treated as a perturbation of So- 

In ordinary perturbation theory, we keep the unperturbed part in the exponential and ex¬ 
pand the perturbation into a series. For convenience, we introduce a formal coupling constant 
g and rewrite the action as S = So + gS i. Then we obtain the ordinary perturbation series 
for the free energy with respect to g as follows: 



( 2 . 2 ) 


We can obtain the free energy up to order n in the perturbation theory by truncating the 
infinite series up to order n and setting g — 1. 

In general, the perturbation series has a finite (or maybe even zero) convergence radius 
with respect to the parameters appearing in the action, for example, the inverse mass squared 


1 /ml in So- To see what happens and when the perturbation theory fails, let us consider a 
zero-dimensional 0 4 one-matrix model. 


2.2 The 0 4 matrix model 


The action of the zero dimensional (j) A matrix model is 


S = + iTY0 4 , 


(2.3) 
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where 0 is an TV x TV hermitian matrix. In the ordinary perturbation theory, we treat the 
quadratic part, 7?t.q0 2 / 2, as the unperturbed action So and the 0 4 /4 part as the perturbation 
Si. Here, as an example, we consider the free energy. After introducing a formal coupling 
constant g and expanding the exponential part with respect to it, we obtain 



and we can estimate this series using the ordinary Feynman diagram technique. Here we 
consider the large -N limit, where only planar graphs contribute to the free energy. Thus, in 
the ordinary perturbation theory, the free energy is given by a series with respect to g/m^ as 
follows: 


f =o 


=-In 

2 



0 ( 9 3 


(2.5) 


On the other hand, one can evaluate this free energy exactly by analyzing the eigenvalue 
distribution |L5|, and it is known that the convergence radius of this series is 1/12. Thus, 


after we set the formal coupling g to 1, this series converges when 7 Uq >12. This means that 
we cannot calculate the free energy for the massless case, i.e. m q = 0, as the limit of this 
perturbation series. Because the IIB matrix model of interest does not have a quadratic term, 
it corresponds to the massless case, in which the ordinary perturbation theory does not work. 


To overcome this difficulty we introduce a new method, the improved mean held approxi¬ 
mation, and obtain an improved perturbation series. 


2.3 Improved mean field approximation 

When we apply the improved mean held approximation scheme, we hrst introduce a “mean 
held” S m (x, a) which can be easily integrated, for instance, a quadratic term. Here, a is a set 
of parameters in the mean held, and we will tune it later to make the approximation better. 
Then, we rewrite the original action as 

S^S 0 + S 1 =* S = S m + (S 0 + S 1 -S m ), (2.6) 

and we take S m as an unperturbed action and the part within parenthesis as the perturbation. 
We introduce a formal coupling constant g as before. Thus the action becomes S = S m +g(So+ 
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Si — S m ), and the expansion of the exponential with respect to g yields another perturbation 
series, which we call the improved perturbation series. Finally, we tune the set of parameters 
a to make the improved series converge, as we explain below. We call this procedure the 
“improved mean field approximation.” In particular, if we consider the first order of this 
approximation for the free energy defined in (|2.5|) , and tune the parameter a using the condition 
dF/da = 0, it is simply the ordinary mean field approximation. 

As an example, we consider again the zero-dimensional 0 4 matrix model. 


2.4 IMFA for the 0 4 matrix model 


In this case, we consider a quadratic term as the mean held: 


S m = — Tv( p 


(2.7) 


Then, we construct a modified action, including a formal coupling g , as 


m 


rrin 


S =—Tr (f) z + g — -Trip 2 + -TrA 


m 


-Tr0" 


= 9 4 + m 2 + g (m 2 -m 2 ) 2 

4 


( 2 . 8 ) 


We emphasize that we do not need to calculate graphs again. By comparing ( |2.8|) with the 
original action Q2.3|) , we have only to substitute m 2 + g(m/ ) — m 2 ) for the mass squared mg 
appearing in the ordinary perturbation theory. 

Therefore, by substituting m 2 + g(ml — m 2 ) for tuq, and re-expanding with respect to g, 


flT5|) becomes 

-improved 


111 


+ 


9 


9 

- o 9 


m 2 + g(ml — m 2 ) J 2 ( m 2 + g[m q — m 2 )) 2 

1 


8 (rn 2 + g(ml — m 2 )Y 


+ 0{g i 


= -In 

2 


m“ ) 2 V nr 4 


i 2 2 

1 771 q — m 


+ 


m “ 
2 


9 1 ml — m 2 1 (ml — m 2 ) 2 


^ 8 


m v 




+ 0(g 3 


(2.9) 


In this form, the massless limit is no longer singular and can be taken simply by setting 
ml = 0. Furthermore, from this form, we know how to calculate the improved perturbation 
series even for a massless theory, to which the ordinary perturbation theory cannot be applied. 
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First, we add a mass term m 2 </> 2 /2 to the original action by hand and calculate the free energy 
perturbatively. Then we substitute m 2 — gm 2 for m 2 and re-expand it with respect to g. 
Finally, by setting g — 1, we obtain the improved perturbation series for the massless theory. 
As we see, this prescription is useful for calculating the IIB matrix model, which also has no 
quadratic term. 

Let us return to the case of a general value of m q. The final step of this prescription is to 
tune the parameter m 2 so that the improved series converges. To elucidate this situation, in 
Fig. B we plot the improved free energies with respect to mg at various orders for m 2 = 2 
and 4. It is clear that the improved series converges quite well in some domain. The position 
of the domain of the convergence depends on the choice of the parameter m 2 , for example, 
around mg = 0 for m 2 = 2 and around m g = 4 for m 2 = 4. Here, we have first plotted the 
improved series with respect to mg for various values of m 2 , and then determined the value of 
m 2 to make the series converge around the value of mg in question. 




Figure 1: Free energies with m 2 = 2 (left) and m 2 = 4 (right). The 
horizontal axis denotes mg. 


2.5 The plateau and how we identify it 

There is an easier way to find an appropriate value of the parameter m 2 than that considered 
above for a fixed value of the bare parameter mg, for example, the massless case mg = 0. Fig. 

displays the improved free energies as functions of m 2 with mg = 0. We see that a plateau 
develops at higher orders. Intuitively, this fact is very natural, because m 2 is an artificially 
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introduced parameter, and the true value of the free energy must be independent of its choice. 
Some years ago, Dhar and Stevenson advocated a “principle of minimal sensitivity” [14]. They 
stated that in the improved perturbation theory one should choose the set of parameters such 
that the improved quantity is stationary with respect to it, i.e. dF imwoveA /dm 1 2 = 0 in our 
case. However, we claim that the criterion for good approximation should be the existence of a 
plateau; that is, one should choose the set of parameters to be on the plateau. This claim was 
first made in a previous paper |T|, and some people have been studying the properties of the 
plateau and seeking a better definition, especially in the case that there are many parameters 
in the mean fields. In such cases, it becomes difficult to identify a plateau, because we cannot 
visualize it easily.^ 






Figure 2: Free energies with = 0 at 7th, 8th, 27th order and from 
Oth to 29th order. The horizontal axis denotes m 2 . 


To obtain a concrete criterion for the plateau, we plot the extrema of the improved free 

1 Nishimura, Okubo and Sugino proposed “the histogram prescription” to identify a plateau for the many 

parameter case jl3| . 
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energy of the massless 0 4 matrix model with respect to m 2 in Fig. It can be seen from this 
figure that there are two characteristic extrema up to 10th order, except at the two lowest 
orders. At each order, one of the two extrema gives the highest free energy, which is quite 
stable, and the other gives the lowest free energy, which decreases as we go to higher orders. 
Hereafter we call the highest and the lowest extremum the “overshooting” and “undershooting” 
extremum, respectively. The other extrema tend to accumulate between these two extrema, 
and they are expected to form a plateau. In fact, if we take a closer look at a plateau, we 
find that it consists of many extrema. In other words, the accumulation of extrema leads to 
a plateau. In general, at much higher orders, a finite number of extrema yield higher or lower 
free energies. However, the situation remains unchanged in which many of the other extrema 
accumulate between these higher and lower extrema and contribute to the formation of the 
plateau. 



Figure 3: Extrema of the improved free energy of the massless </> 4 
matrix model from 1st to 10th order. 


This is the key observation for finding a plateau presented in this paper. In the next 
section, we explore the extrema of the free energy of the IIB matrix model and see whether 
they accumulate to form a plateau or not. 



3 Free Energies of IIB Matrix Model 


In this section, we would like to analyze the IIB matrix model [[|. The action is given by 

S„b= —Tr( (3-1) 

where A M {/i — 1, • • • , 10) and ip 01 (a = 1, • • • , 16) are all N x N hermitian matrices trans¬ 
forming as the vector and left-handed spinor representations under SO(IO). go is the only 
parameter in this model, which is dimensional, and we take the large N limit while fixing 
the ’t Hooft coupling g^N to 1. This is defined as a zero-dimensional reduced model of a 
ten-dimensional super Yang-Mills theory. Because it has no quadratic term, we cannot apply 
the ordinary perturbation theory. 

We use the improved mean field approximation to evaluate the free energy. Following the 
method explained in the last section, we add and subtract a quadratic term for bosonic and 
fermionic matrices as mean fields. Then our action becomes 

S = - N TV ( | [A„ AJ 2 - A i }T“lA m <!>])+S m - gSn, (3.2) 

Sm =\ (C-^N^A^) - (Cf l ) m N TV (V“<). (3-3) 

where S m is the quadratic term introduced as the mean field^]. Here we introduce the formal 
coupling constant g as g\ in (|3.1|) . The coefficients C iJiV and^ = u IJil/p Y ,wp /3\ are the propagators 
for bosonic and fermionic matrices in the perturbation theory, respectively. C^ v ) is a second 
rank symmetric tensor, and U[ t w P \ is a third rank antisymmetric tensor. 

Here we should comment on how these mean fields are constructed. The symmetries of the 
original IIB matrix model are the matrix rotation U(IV), ten-dimensional Lorentz symmetry 
SO(IO), translational symmetry, which changes A^ to A^ + const, x 1^, and the type IIB 
supersymmetry. Our mean field preserves U(iV), while it breaks SO(IO) and supersymmetry. 
One might worry that this leads to an inconsistency because the existence of the type IIB 
supersymmetry plays an important role in the IIB matrix model. However, if the supersym¬ 
metry is restored in the true ground state of our model, the parameters will go to C = oo 
and u = oo and our mean field should vanish. This is the standard story of models in which 
dynamical symmetry breakdown occurs, as in the Nambu-Jona-Lasinio model. In this case, 
even if we introduce a mass term that breaks the chiral symmetry, it vanishes in the phase 
where the symmetry is restored. If C and u are still finite after the large N limit is taken, 

2 C is the charge conjugation matrix defined as t C = —C, CT M = — Y (g = 1, • • • , 10). 
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both the Lorentz symmetry and supersymmetry are spontaneously broken. This is the very 
scenario we expect. 

Here we summarize the procedure we carry out below. Following the prescription for the 
massless case described in § |2.4| , we first calculate the free energy perturbatively for the action 

S' = - iVTr( | {A„,A„} 2 - Ajr ) + S m . (3.4) 

Then we replace C _1 and u l with (1 — g)C~ 1 and (1 — g)u~ l , respectively, in order to 
recover the contributions from the —gS m term. Finally, we obtain the improved free energy 
by expanding the result with respect to g. 


3.1 The 2PI free energy and ansatz 

To obtain the free energy, we only have to calculate planar connected vacuum graphs. However, 
in such calculations, we have to treat too many graphs at higher orders. In order to avoid 
this problem, we introduce a two-particle irreducible (2PI) free energy. The 2PI free energy is 
considered in Ref. Jl[| and it is deeply related to the Schwinger-Dyson equations, as discussed 
in Ref. [|lj. In this paper, however, we use the 2PI free energy only as a tool to obtain the 
ordinary free energy easily. From this viewpoint, the definition and properties of the 2PI free 
energy are summarized as follows: 

• The 2PI graph is the graph that contains no self-energy graphs as its subgraphs; that 
is, it is two-particle irreducible. This means that propagators in a 2PI graph can be 
regarded as the exact propagators. 

• Suppose G is a sum of planar, vacuum, connected 2PI graphs in some theory. Then the 
ordinary free energy in the planar limit is given by the Legendre transformation of G 
with respect to the exact propagators. 


In the IIB matrix model, we also force G to contain no tadpole graphs, because all the one- 
point functions satisfy ( A M } = 0, due to Lorentz invariance. 


Once we obtain a 2PI free energy G of (|3.4j) up to some order, the next task is to perform a 
Legendre transformation. At this stage, however, we face the new problem that there are too 
many parameters to carry out the transformation. Therefore we need some ansatz that reduces 
the number of parameters to a tractable one. We use the same ansatz as in Ref. [jfj; that is, 
we assume several unbroken Lorentz symmetries that restrict the mean held parameters 
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and u^p. Two peculiar examples, which we consider in this paper, are the SO(7) x SO(3) 
case and the SO(4) x SO(3) x SO(3) x Z 2 case. As we comment below, these two ansatz are 
of particular importance among other possibilities. When the unbroken symmetry is SO(7) 
x SO(3), which we call the SO(7) ansatz, the parameters associated with the bosonic fields 
Ap are limited as C fW = diag(Vi, Vi, Vi, V t , Vi, Vi, Vi, V 2 , V 2 , V 2 ), and those with the fermionic 
fields il) a are restricted according to Ms, 9,10 = —m 9 ,s,io = (cyclic) = u, while the others are 
zero. The other one is called the SO(4) ansatz, which preserves SO(4) x SO (3) x SO(3) x 
Z 2 symmetry. The SO(4) symmetry acts on the indices /i = 1, ,4, and the two SO(3) 

symmetries act on g = 5, 6, 7 and g = 8, 9,10, respectively. The Z 2 symmetry exchanges 
g = 5, 6, 7 and g — 8, 9,10 directions, that reverses the 1st direction so that it is an element 
of the SO(IO). 

Then, in the SO(4) ansatz, the parameters arc limited as 
C, w = diag(Vi, Vi, Vi, Vi, V 2 , V 2 , V 2 , V 2 , V 2 , V 2 ), m 5 , 6 , 7 = m 8 , 9 ,io = u/\/2, up to the cyclic permu¬ 
tation of indices with signature, while the other components are zero. 

Here we should comment on the other ansatz. In Ref. |L|, we consider, besides the ansatz 
mentioned above, SO(l), SO(2), SO(3), SO(5) and SO(6) ansatz. For each ansatz, SO(n) 
represents the residual symmetry for the expanded direction, which at the end should be 
understood as the space-time dimension. According to the analysis given in Ref. [;l|, the 
SO(2) and SO(3) ansatz restore SO(4) symmetry, the SO(5) and SO(6) ansatz restore SO(7) 
symmetry, and the SO(l) ansatz has no extrema and, of course, no plateau. (See Ref. flU 
for details.) After all, we concentrate on the SO(4) and SO(7) cases, which behave quite 
differently, as we now see. 


3.2 Calculation of the free energy and the extent of space-time 


Now we have all the tools needed to compute the improved free energy. Our first task is to 
calculate the 2PI free energy G for the action (|3.4|). It can be computed as 


G 

N* 



( 1 


© 


+ .9 2 i; 



- ^ ln(det C) + ^ln(det^) + g ^-^(tr^C' 2 ) - (tr^C) 2 ) - ^Cp U tr a (^T^T 
+ g 2 ^(tr M (C 4 ) - (tr,d 2 ) 2 ) + i^C pA tr Q (^r^r^r A )] + (D(g 3 ), 


(3.5) 
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where the solid and dashed lines represent fermion and boson propagators, respectively, tr p 
represents the trace in the vector representation of SO(IO), and tr Q is the trace taken in the 
left-handed spinor representation of SO(IO). In order to count the order, the formal coupling 
constant g is explicitly inserted. We calculated this 2PI free energy to 7th order in g. The 
numbers of graphs we computed are 2 at Oth, 1st and 2nd orders, 4 at 3rd order, 12 at the 4th 
order and 49 at 5th order. All of these were calculated in Ref. []T[ and all the explicit graphs 
are presented in that paper. Our new result is for 6th and 7th order. The numbers of the 
graphs at 6th and 7th orders are 321 and 2346, respectively. The generation and calculation 
of the planar 2PI graphs have been now totally automatized^. 

Next, we perform the Legendre transformation. We define the variables conjugate to Vi 
and u as 

M>= wV v,u) ’ (3 ' 6) 

m =-^-G(V,u), (3.7) 

and obtain the ordinary free energy via the Legendre transformation as follows: 

F{M\m) = G -J2 Mi Vi~mu. (3.8) 


Finally, the improved free energy Fi mpr oved to order g k can be obtained by subtraction, as 
in the case of the massless q i 3 4 model discussed in (Q: 

^improved = ~ 9 M *> m ~ 9 m ) \k • ( 3 ‘ 9 ) 

Here |& denotes expansion with respect to g, ignoring 0(g k+1 ) terms and setting g — 1. 

Thus we obtain the improved free energies to 7th order. As we have discussed, in general 
it is very difficult to identify where plateaus are and how they grow. 

In the next subsection, we find the extrema of these free energies at several orders, and we 

compare their distribution with that for the matrix model. 

3 In order to check the algorithm for generating graphs, we have used a matrix b 4 -QED model that provides 

the same 2PI graphs as those of the IIB matrix model and is easier to compute analytically. Then, we confirmed 

that the set of graphs and the symmetry factors are all correct. (See Appendix |A|.) 
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3.3 Extrema of the improved free energies and the extent of space- 
time 

We now have the improved free energies of the SO(4) and SO(7) ansatz to 7th order in the 
improved perturbation. In order to search for a plateau, we list all the extrema of these free 
energies, as done in The extrema we have found are listed in Table |l|. 

Now we consider the “extent” of space-time in two directions, which is defined by the 
moment of the eigenvalue distribution as follows: 

1 BF 

« 2 = ( - Tv A, 2 ) = ~g UT (M i - gM\m- gm) . (3.10) 

1 BF 

r 2 = ( - TvA m 2 ) = - gM\ m-gm). (3.11) 

We call R the extent of “our” space-time and r that of the internal one. The values of R, r 
and their ratio, p = R/r , are also shown in Table |l] for each extremum of the SO(4) and SO(7) 
ansatz. 

In Figs. [I] and |5|. these extrema and the ratios p are plotted for each ansatz. 

We find that these two ansatz yield significantly different behavior. The SO(7) ansatz has 
fewer extrema than the SO(4) ansatz, and, in particular, it has no extrema at even orders. 
Therefore, we speculate that extrema for the SO(7) ansatz do not accumulate to form a 
plateau. This means that the SO(7) ansatz is not a realistic assumption for the eigenvalue 
distribution of the IIB matrix model, and thus a seven-dimensional flat space-time would not 
be realized as its stable vacuum. Even if the SO(7) ansatz does develop a plateau at higher 
order, it would not be regarded as a compactification, because the ratio of the extents seems 
to be stabilized around 2. 

In contrast, the SO(4) ansatz exhibits interesting behavior. In this case, the number of 
extrema grows as the order increases. It seems that even orders are not stable compared to 
odd orders. We observe similar situations in various models. Actually, in the zero-dimensional 
cj) A model, even lower orders have no extrema, while odd orders develop a plateau even at lower 
orders |IJ . Assuming this is the case, we ignore even orders. Then we find peculiar behavior: 
There are two characteristic extrema that can be taken as the counterparts of the overshooting 
and undershooting extrema of the 0 4 matrix model, and they are unrelated to a plateau. It 
seems that the other extrema tend to accumulate and can be expected to form a plateau. We 
further observe that the ratio of the extents for the SO(4) ansatz are around 2 or 3 for these 
undershooting and overshooting extrema, whereas the other extrema have rather large ratios, 
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ansatz 

order 

F 

p = R/r 

R 2 

r 2 

S0(7) 

1st 

5.52272 

1.95530 

0.410551 

0.107383 

3rd 

5.62672 

1.91133 

0.463883 

0.126981 

1.62094 

2.15683 

0.588327 

0.12647 

5th 

5.52146 

1.92523 

0.481755 

0.129975 

3.00676 

2.19364 

0.573357 

0.11915 

1.49243 

2.0943 

0.648315 

0.147812 

7th 

5.45127 

1.93442 

0.491098 

0.13124 

3.89291 

2.39213 

0.640411 

0.111916 

1.92312 

1.62093 

0.33974 

0.129306 

1.93939 

1.93483 

0.519378 

0.138739 

SO (4) 

1st 

6.1533 

1.85728 

0.562580 

0.163090 

3rd 

6.34486 

1.85336 

0.650887 

0.189489 

0.696885 

3.05943 

1.34914 

0.144137 

4th 

1.17141 

4.42435 

1.90023 

0.0970747 

1.55668 

2.58548 

1.1568 

0.173051 

5th 

6.26112 

1.92138 

0.700007 

0.189617 

2.85746 

5.00919 

2.12045 

0.0845069 

0.5426 

3.75008 

1.81111 

0.128784 

0.0776919 

6.34998 

2.92093 

0.0724395 

0.473299 

3.31357 

1.65917 

0.151112 

6th 

-2.62258 

5.45278 

2.69232 

0.0905505 

-2.78918 

7.34015 

3.55225 

0.0659317 

-0.642089 

7.36892 

3.05248 

0.056214 

0.0100066 

3.50285 

1.30668 

0.106494 

7th 

-2.66095 

3.11512 

1.17719 

0.377896 

6.19216 

1.97197 

0.73242 

0.188347 

5.16286 

6.42722 

2.84737 

0.0689285 

2.75812 

6.8741 

2.22076 

0.0469969 

3.77441 

9.26784 

3.82752 

0.0445615 


Table 1: Extrema of the free energy and the extent of space-time for the S0(4) and S0(7) 
ansatz. 
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Figure 4: Extrema of the free energy (upper) and ratio of the extents 
(lower) for the SO(7) ansatz. 
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Figure 5: Extrema of the free energy (upper) and ratio of the extents 
(lower) for the SO(4) ansatz. The numbers assigned at 7th 
order display a correspondence between the extrema and 
the value of the ratio of the extents. 
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around 6-9. This implies that on the plateau, the ratio of the extents takes a large value. In 
short, we can conclude that the SO(4) ansatz develops a plateau, and it predicts a quite large 
value for the ratio of the sizes of the internal and external directions. This indicates that onr 
scenario for spontaneous compactification to a flat four-dimensional space-time is promising. 

4 Conclusions and Discussions 

We have performed the improved mean field approximation for the IIB matrix model up to 
7th order and obtained the following conclusions: 

• We first conclude that the SO(7) and SO(4) ansatz exhibit different types of behavior, 
as described below. 

• The SO (7) ansatz has fewer extrema than the SO(4) ansatz, and it does not have a 
tendency to form a plateau. The eigenvalue distribution of the SO(7) ansatz is rather 
isotropic, and will not be realized as a compactification vacuum even if it has a plateau. 

• The SO (4) ansatz has many extrema at higher orders and, except for the two special ex¬ 
trema mentioned below, it yields a quite large ratio of the extents of the four-dimensional 
space-time and that of the internal one. 

• The extrema of the SO(4) ansatz are distributed, as in the 0 4 matrix model, which 
develops a plateau. The 0 4 matrix model has an overshooting and an undershooting 
extremum, which are located over and under the plateau, and the other extrema tend 
to accumulate. Thus, we conjecture that the IIB matrix model exhibits similar behavior 
and develops a plateau under the SO(4) ansatz. 

At this stage we do not have clear plateaus for any ansatz, and therefore we cannot defi¬ 
nitely tell which vacuum is realized in the IIB matrix model. Indeed, in order to do this, we 
need to identify the plateau for each ansatz, if such exists, and compare the values of the free 
energies at the plateaus. If an ansatz has no plateaus, we conclude that it is not realized as 
a vacuum. We expect that our SO(4) ansatz is close to the true vacuum, which reproduces 
our universe, and it has a plateau where the free energy has the lowest value. In order to con¬ 
firm this, we should analyze higher orders in the improved perturbation series. As mentioned 
above, our calculation is now totally automatized, and it seems possible to carry out further 
analysis with the help of supercomputers. 
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A The 2PI free energy of </> 4 -QED type matrix model 


In this appendix we consider a matrix model that has a 0 4 type and a QED type interaction 
term. The action is given by 


S = JVTr | #A 2 + ^A 1 + JA + #V> 2 + gAip 2 


(A.l) 


where A and # are N x N matrices and # is assumed to have flavor /. As we discuss below, 
we have introduced a source J for A in order to cancel the tadpole graphs. 

By comparing this model with the IIB matrix model with the mean held action ( |3.4|) , it 
is easy to see that ( |A.1| ) generates the same vacuum graphs with the same symmetry factors 
as the IIB matrix model. These two models have the same types of propagators and vertices, 
except for the source term JA, which we explain below. Here mg and mj? play the roles of 
C~ x and ^ _1 , respectively, and / should be set to —1, because # is fermionic in the IIB matrix 
model. However, we do not fix /, in order to classify graphs via the number of #-loops. It 
is worth noting that in (|A.1|) there is no symmetry that forces the one-point function ( A ) 
to vanish. This is in contrast to the situation in the IIB matrix model, where it cannot have 
a non-zero value, due to the Lorentz symmetry. In order to eliminate unwanted one-point 
functions from our model, we choose J in such a way that ( TrA ) = 0 order by order in the 
perturbation theory with respect to g. 

Our aim is to confirm that our list of the planar 2PI graphs of the IIB matrix model 
is complete. For each graph, we read off the number of boson propagators (B), fermion 
propagators (F), fermion loops (L), A 4 vertices (V) and A# 2 vertices (Y), and we deduce the 
2PI free energy of (|A.1|) as 

1 


symmetry factor 




Y-B-F f L 


m. 


7 , 


(A.2) 
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where the summation is taken over all possible planar 2PI graphs without tadpoles up to the 
given order of g. Note that in this expression cancellation between different graphs never 
occurs, because the symmetry factor is always positive. We compare this with the 2P1 free 
energy of ( |A.1| ) in the large -N limit, computed with a completely different method that we 
now explain. 


A.l Loop equation 


The field redefinition in (|A.1|) gives 


S 



(A.3) 


In order to recover the original parameters in (|A.1|) , we only need to do is to make the 
replacements g —> g/m B , J —> J/y/m B , A —> yjm B /m B ■ 

For the purpose of computing the 2P1 free energy of (|A.3|) , we first compute the two-point 
function ( Tr(A 2 ) )/N in the large-IV limit, retrieve the parameters m B and m?, and integrate 
it with respect to to obtain the ordinary free energy. Finally, we carry out its Legendre 
transformation in terms of m B and mg to obtain the 2PI free energy. 


Performing the Gaussian integration in ( A.3 ) with respect to ift, we obtain 

/ ^k^k (-A g) p+q 


S, s = NTr (A 2 + La 4 + JA) - £££ 


2 p + q 

p = 1 9=0 ^ 1 


- p+q C q Tb(kP)Tb(kP) 


(A.4) 


In order to compute the two-point function ( Tr(A 2 ) )/N of this model, we start with the loop 
equation (Schwinger-Dyson equation) 

°= f dA-^Tr(A n t a )e- s ^, (A.5) 


where t a is the orthogonal basis of hermitian N x N matrices: Tr= 5 a/3 . By defining 


w n =-(Tr(A n )), 


(A.6) 


this equation gives us the relation between these correlation functions, 

n —2 oo oo 

W n ^ ^ '^rn'^n—2—m Q ^n+2 J^n— 1 “1“ f y^(-A S rv,- \C q Wn+p—2'^qi (-A.7) 

m =0 p =1 9=0 

where we have used the factorization property in the large -N limit. Note that this equation 
holds even in the n — 1 case if we ignore the first term on the right-hand side. The loop 
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equation enables us to determine w n order by order in the perturbation theory as follows. 
First we expand each w n in terms of g : 


VJn = 

k =0 


(A.8) 


Substituting this into 
following equation: 


(A.7) and comparing both sides to each order in g, we obtain the 


n —2 k 

w ( k ) — \ \ _ (A: 2) _ t (k) 

w n ~ /_^ w m W n-2-m w n+2 J w n- 1 

m= 0 Z=0 

k k—p k—p—q 

+ /EE (A.9) 

p= 1 <j=0 z=o 


Because w o = l/iV( Trl ) = 1, we have the “boundary condition" 


Wn 


= 


Note that 


only the quantities Wm with l smaller than or equal to k appear on the right-hand side and 
that when Wm' 1 appears on the right-hand side, m is always smaller than n. Because of this 
property, (|A.9|) together with the boundary condition determines all the . 


Once we determine vj\ k> up to a given order, for example, the 14th, corresponding to the 

(k) 

7th in the 11B matrix model, we can tune J order by order in such a way that w\ vanishes 
for each k. Substituting the resulting expression for J into w 2 , we can derive the two-point 
function up to the order we desire. Because of this procedure, graphs containing a tadpole as 
a subgraph no longer contribute to w 2 k \ and therefore the free energy obtained by integrating 
w 2 = Thk9 kw 2 ^ gives the sum of all the planar vacuum graphs without tadpoles. Once we 
have the free energy, it is easy to carry out the Legendre transformation to obtain the 2PI free 
energy. 


We have compared the two forms of the 2PI free energy obtained using the two totally 
different methods and found complete agreement. This proves that our list of all planar 2PI 
graphs of the 11B matrix model is complete. 
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